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ABSTRACT 

We present new calculations of unified line profiles for hydrogen perturbed by 
collisions with protons. We report on new calculations of the potential energies and 
dipole moments which allow the evaluation of profiles for the lines of the Lyman series 
up to Lyman 8 and the Balmer series up to Balmer 10. Unified calculations only 
existed for the lines Lyman a to Lyman 7 and Balmer a including the quasi¬ 
molecule. These data are available as online material accompanying this paper and 
should be included in atmosphere models, in place of the Stark effect of protons, since 
the quasi-molecular contributions cause not only satellites, but large asymmetries that 
are unaccounted for in models that assume Stark broadening of electrons and protons 
are equal. 
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1 INTRODUCTION 

Collisional broadening occurs whenever a radiating parti¬ 
cle interacts with a perturber during the process of emis¬ 
sion or absorption of radiation. This interaction is known 
to cause not only a broadening of the line profile, but also 
the possibility of appearance of satellite-lines. Satellites are 
explained by the formation, during the collision, of a tran¬ 
sient bond state which has energy levels other than those 
of the non-interacting particle, producing a spectrum unlike 
the obtained for the isolated particle. 

Hydrogen satellite-lines were already observed in some 
lines of the Lyman series in stars, and in the Lyman and the 
Balmer series in laboratory plasmas. In stars, these detec¬ 
tions are mainly in spectra of white dwarf stars, but are also 
reported in A-type Horizontal Branch stars (which may be 
extremely low mass white dwarfs, ELMs) and A Bootis stars 
l|Allard et al.ll200cil b The first detections were two satellites 
in the red win g of Lyman a, close to 1400 (lGreensteinlll980l : 
IWeenerl li 9821 ) and 1600 A (iHolm et al .11 1 9851) . of DA-type 
white dwarfs with effective temperatures around 15 000 K. 
As they could not be explained by transitions of the hydro¬ 
gen atom, which are the only lines expected to b e observed 
in DA whi te dwarfs, they w e re exp lained bv lKoester et all 
Jl985h and iNelan fe Wegner! (Il985l l as quasi-molecular con¬ 
tributions due to H^ and H2, respectively. 

Since then, the effect of said two contributions has 
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been included in atmosphere mode ls, offering a good diag¬ 
nostics of stellar parameters (e. g. lAllard fe Koesten Il992l : 
lAllard. Kielkopf fe Loeilletl 120041 ') . With the improvement 
of observational techniques and increase in data qual¬ 
ity, sa tellites in othe r spec tral lines have also been iden¬ 
tified. iKoester et alJ Jl996h reported two strong absorp¬ 
tions in the red wing of Lyman /3, at 1060 and 1078 A, 
in the spectra of the DA white dwarf Wolf 1346 ob¬ 
tained with the Hopkins Ultraviolet Telescope (HUT). In 
the same paper they already presented a model for the 
absorptions assuming a constant dipol e mom ent for the 
transitions. lAllard. Kielkopf fe Feautriei! Jl998h approached 
the problem again, now consideri ng the true variat ion of 
the dipole moment. In Lyman 7, IWolff et alJ (1200il l iden¬ 
tified a satellite with wavelength around 990 A in the 
Far Ultraviolet Spectroscopic Explorer (FUSE) spectrum 
of the star CD -38° 10980. Models were evaluated by 
Allard et al. <12004 1. assum i ng co nstant dipole moments, and 


All ard. Noselid ze & Kruk d2009f) . considering their depen¬ 


dence with distance. 


In laboratory plasmas, experiments focused on quasi- 
molecular contributions were done for Lyman a and the 
Balme r series. Lyman a was studied by iKielkopf fe Allardl 
dl995l) using laser-produced plasma in pure hydrogen. They 
detected five satellites: the already known features at 1400 
and 1600 A and three others closer to the line centre, at 
1234, 1172 and 1268 A, explained by their models as con¬ 
tributions due to the quasi-molecules Hj (1234 A) and H2 
(1172 and 1268 A), but not detected in white dwarfs spec- 
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tra because the l arge pressure b roadening of the sp ectral 


l arge pressure b roadening ot t ne sp e 
lines hides them. iKielkopf. Allard fc Decrettel d2002T) stud¬ 
ied Balmer a and identihed six satellites in the theoretical 
wing, of which only a stronger absorption at 8650 A is ob¬ 
served in the plasma spectrum. For Balmer p and Balmer 7 , 
on the other hand, no satellites were identified, because only 
the line centre exp e rimen tal profiles were published so far 
llFalcon et al.l 120121 . 12014 ) and usually there are no satel¬ 
lites in this region, which corresponds to radiation_when_the 
particle s are far apart. In these papers and in iFalcon et al.l 
d 2010 ah . they use x-rays generated from a z-pinch dynamic 
hohlraum to drive plasma formation in a macroscopic gas 
cell and generate absorption spectra for the Balmer lines, 
with the intention of benchmarking white dwarf stars at¬ 
mosphere models, though the published data are scarce to 
date. 

Accuracy in model atmospheres for white dwarf stars is 
essential for a good determination of their physical parame¬ 
ters, especially effective temperature and surface gravity, as 
the fit of synthetic spectra is t he most successfu l method to 
obtain those parameters (e.g. iTremblav fc Bereeronl f2009h . 
White dwarf stars are the final evolutionary state of all 
stars with masses be l ow 8-1 0. 5 Mq , dependin g on metallic- 


_e pe 
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120081 : 1 Doherty et al. 20l4 , corresponding to over 97 % of 
the total number of stars. They generally do not present 
ongoing nuclear fusion and only emit their residual ther¬ 
mal energy, resulting in a slow cooling of the star. Because 
they have radii of the same order of the Earth’s radius, 
they have a small surface area, so their cooling rate is small 
and therefore their cooling times are large; it takes approx¬ 
imately 10 10 years for the effective temperature of a nor¬ 
mal mass white dwarf to decrease from 100 000 K to near 
5 000 K. As a result, the cool ones are among the oldest ob¬ 
jects in the Galaxy, being extremely useful in the study of 
stella r formation and ev o lution history of the Milky Way 


(e.g. \\ ingel et al. 19871: Be n. reron, ^Jaffer^JJebert 


Liebert, Bergeron fc Holberd 120051 : iMoehler fc Bono 
Tremblay et al.ll20l4 . 


1992; 


2008; 


Currently the atmosphere models in use to fit white 
dwarf stars and estimate their physical parameters only 
take into account quasi-molecular contributions for a few 
lines of the Lyman series, when they consider such con- 
tributions at all: Ly man a, Lv man p and Lyman 7 (e.g. 
lAllard. Noselidze fcKrufl 120091 ). For the other lines of the 
Lyman series and the Balmer series, the interaction between 
hydrogen atom and free protons is approximated simply as 
electron Stark effect, which causes only a broadening that 
adds up to the effect of the electrons. The quasi-molecular 
effects are thus neglected, disregarding that they are respon¬ 
sible not only for a number of known satellites but also for 
strong asymmetries between the red and blue wings of all 
lines. The far red wing of Lyman a (A > 2 500 A), for exam¬ 
ple, turned out to be the missing opacity in the blue part of 
the optical spectrum of cool white dwarfs (Teft < 6 000 K) 
of pure hydrogen composition. The inclusion of the broad¬ 
ening by collisions with H 2 was the key ele men t in obtaining 
such result, as shown by I Kowalski fc Saumonl (|2006j ). How¬ 
ever, contributions to other lines remain unaccounted for, 
despite the fact that their inclusion may be important to 
obtain accurate values for effective temperature and surface 
gravity. 


This lack of accuracy results, for example, in a discrep¬ 
ancy in the calculated mean masses of DA white dwarf s 
when distinct methods are applied. IFalcon et al.l (feoiQbl ). 
analysing 449 DA white dwarfs from the European South¬ 
ern Observatory Supernova la Progenitor Survey (SPY), via 
the gravitational redshift method, obtained a mean mass of 
(Mda) = 0.647lt§;oi4 Mq. A consistent v alue of (M p A ) == 
0.630 ± 0.028 Mq was obtained by iRomero et al.l (|2012h 
with seismology of 44 ZZ Ceti white dwarfs. On the other 
hand, when analysing 2 217 DA stars from the Sloan Dig¬ 
ital Sky Survey (SDSS) data release 7 with signal-to-noise 
ratio larger than 15 and T e g > 13 000 K [to avoi d the sys¬ 
tematic super-e st imati ve reported bv lKepler et alJ |2006)], 
iKleinman et al.l (120131 ) found a mean mass of (M D a) = 
0.593 ± 0.002 Mq, which is sig n ificantl y sm aller than the 
values obtained bv IFalcon et al. (I2010bl4 and IRomero et all 
ll2012ll . Applying correct i ons to 3D convection , as c alcu- 
lated by ITremblav et al.l (l2013l h iKepler et al.l (120151 ') ob¬ 
tained a mean mass of 0.662 ± 0.003 Mq in the analy¬ 
sis of 1 659 DA white dwarfs from the SDSS data release 

10 with signal-to-noise ratio larger than 15 and T e g_> 

10 000 K , whi ch is a value closer to the results of lFalcon et alJ 
(l 2010 bh and IRomero et alJ (120121) . Hence the inconsisten¬ 
cies between the spectroscopic method and other meth¬ 
ods are reduced when corrections to 3D convection are ap¬ 
plied, b ut are st ill with out a comp lete solution at present 
dGene st-Beaulicu & Bcr geronll20l4 . 

Falcon et al.l |2014) noted that the inferred electronic 
density of a plasma when fitting H 7 was ~ 40 % smaller 
than the density obtained when fitting H /3, regardless of 
the set of models used. They suggested that this inconsis¬ 
tency would lead to lower determinations of surface grav¬ 
ities and, by consequence, a lower mean mass would be 
obtained. The probable cause for the different values ob¬ 
tained with each line is some missing physics in current 
atmosphere models. In this work we calculate line profiles 
for the quasi-molecular contributions of collisions between 
hydrogen and free protons, which is one of the effects not 
yet consistently included in the models. The potentials and 
dipole mom ents for Hj~ used here are extensions of those 
published by lSantos fc Kepled (120121 ). Our results show that 
the quasi-molecular effects should be included in the atmo¬ 
spheric models, regardless of the detection or not of satel¬ 
lite lines, since they cause also large asymmetries in every 
line. As they are not contemplated in the current calculated 
spectra, they are probably affecting the determinations of 
temperature and gravity for the analysed data. 


2 THEORETICAL BASIS 

There are basically four possible approaches for the calcu¬ 
lation of collisional profiles: the impact theory, the quasi¬ 
static approximation (statistical theory), the unified the¬ 
ory an d the quantum method. The impact theory (e.g. 
iLindholml Il945l ) assumes that the duration At of a colli¬ 
sion is much smaller than the interval between successive 
collisions (~ 2n/Aui, where Au is the angular frequency 
measured from the line centre), being useful for determin¬ 
ing the behaviour of the line centre, but not being able to 
account for the satellites. On the other extreme, there is the 
quasi-static approximation, which assumes At 2n/ Auj, 
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being then a Rood method for obtaining; the line wings (e.g . 
Rohrman. Althaus fe Kepleij l 201 ll : ISantos fe Keplenl 2012 l) . 
Andersonl 1 19521 ') showed that both approaches could be 


obtained from the same formulation, taken the appropri¬ 
ate limits, proposing the unified theory , whic h can give 
the shape of the whole line. iAUard et al„ ( 1199! ?;) developed 
such theory in a quantum mechanical approach, obtain¬ 
ing from it better results for the shape and amplitude of 
the satellites. The pure quantum mechanical method can 
also give the whole line, but, for its implementation, one 
must know all of the states involved in a transition, in¬ 
cluding rovib r ational, deman ding larger computational time 
llZygelman fe Dalgarnol 199C ). We chose the unified theory 
as proposed by Allard et al. ( 19991 ') to the evaluations in this 
work. 

In the unified theory, the normalized line profile mea¬ 
sured from the line centre is given by the Fourier transform 
of an autocorrelation function <F(s): 


I(Alj) = -Re 


$(s)e lA “ s ds. 


(1) 


The autocorrelation function measures the influence of the 
presence of perturber particles - in this case, free protons - 
in the radiatio n of the particle - here, hydrogen. As shown 
by 1 Allard et al.l (|l999i l. it can be written as: 




_ „™ 9 ( s ) 


( 2 ) 


where n is the perturber density and g(s), for a transition 
a = ( i, /) between the initial state i and the final one /, is 
given by 


'[R( o)] 


1 _ n oo noo 

9 <*(s) = ™/ da: 0 / 2nbdbd e 

2^e,e f I ee> > e e / "—oo J 0 

X {e* f ° V '«' lR(t)]dt dU [R(s)] - dee' [-R(O)] } , (3) 

where e and e! denote the energy subspaces containing the 
states which approach, respectively, the initial and final 
atomic states of the transition as R — > oo, with R the dis¬ 
tance between radiator and perturber, d ee / the dipole mo¬ 
ments for the transition, d ee / the dipoles modulated by the 
the Boltzmann factor e _/3Ve , xo the perturber initial posi¬ 
tion, b its impact parameter and V ee i(R) = 14 (R) — V e /(R) 
the potential energy difference between states i and /. 

We adopted this approach to evaluate the line pro¬ 
files I (Aw) for the Lyman series up to Lyman <S and 
to the Balmer series up to Balmer 10. The input data 
for potentials a nd dipole moments were calculated by 
ISantos fe Kenlerl d 2012 lf and are available at their web 
pagfl except for the potential energy data used for the 
evaluation of the line centres, which were obtained with 
an improvement of the code, as will be described in 
Section [ 6 ] Calculations via un ified theory existed in the 
litera ture only fo r Lyman a ( Allard, Kiclkopf jWLoeilletl 
120041 '). L yman B (I Allard, Kielkopf R F eautrieij 19981 '). Ly- 
man 7 ( Allard, Noselidze fe Kr 1 iU~ 20091 ) and Balmer a 
llKielkopf. Allard fe Decrettel l20 o'^L For the oth er lines, cal¬ 
culations were done only by Santos fe Kepled ||2012|) . but 
via the quasi-static approach. Evaluations covering also the 
line centre were in order, not only to improve white dwarf 


1 http://astro.if. ufrgs.br/marcios/index.htm 



Figure 1 . Opacity profiles for every computed line in the Lyman 
series for a temperature of 10 000 K and perturber density of 
10 17 cm -3 . Strong satellites can be identified at the wavelengths 
of 1400, 1080, 1060, 990 and 965 A. The last one has not been 
identified in observations yet. 


atmosphere models but especially to allow co mparison with 
new experim e ntal d ata like those published bv lFalcon et al.l 
ll2010al . T2012I . I 20 I 4 I ). That is the scope of the present work. 


3 LINE PROFILES FOR THE LYMAN SERIES 

The Lyman series corresponds to transitions which have n = 
1 as the final level, for absorption lines, or as the initial level, 
in case of emission. Its wavelengths are in the ultraviolet 
region, between 912 A (Ly 00 ) and 1216 A (Ly a). Here 
we evaluated Ly a, Ly /3, Ly 7 and Ly 5, whose opacity 
profiles are shown in Fig.[T|for a temperature of 10 000 K and 
perturber density of 10 17 cm~ 3 . The behaviour is linearly 
extrapolated for log(<JA) < —26. Many satellites can already 
be identified, as listed in Table [T] but each profiled is shown 
in more detail in the next subsections, for clarity. 


3.1 Lyman a 

Lyman a has contributions due to six transitions of the 
Hj quasi-molecule. The total profile accounting for each 
contribution is shown in Fig. [5] for temperatures of 10 000, 
15 000 and 20 000 K and perturber densities of 10 16 , 10 17 
and 10 18 cnL 3 . Four satellites can be identified, with ap¬ 
proximate wavelengths of 980, 1100 , 1240 and 1400 A. The 
first two are probably too faint to be detected in stellar spec¬ 
tra, at least at these temperatures. The one at 980 A would 
also be blended with the Lyman 7 line in stellar spectra. 
The other two were already observed, as cited in the intro¬ 
duction, in laboratory plasmas and in white dwarfs. We can 
also note that, in large scale, the temperature does not seem 
to affect the profiles, but, in the detailed box, it can be seen 
that the satellite at ~ 1400 A has a consistent behaviour: it 
broadens with the increase in temperature, while its peak is 
displaced to smaller wavelengths. Such behaviour might be 
useful as a good temperature diagnostic for plasmas. The 
effect of increasing density is simply a linear shift (in log 
scale) and a larger broadening. 

In Figs. [3] and [4] we show comparison between observed 
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Table 1. Wavelengths for the maximum of the stronger satellites in each calculated line in the Lyman series 
(for T = 10 000 K and n = 10 17 cm -3 ), the transition(s) that originates them and notes on observation and 
theoretical models. 


Line 

Satellites (A) 

Transition(s) 

Observed? 

Previous model 

Ly a 

979.55 

2p<T„-3d7Tg 

no 

Allard et al. f1994 j* 3, 

Allard et al. l'1994'M 


1118.42 

2sa g -2sa g 

no 


1230.88 

lscr g -4fcr tl 

lscr g -2p7r„ 

ves (Kielkopf &; Allard 1995) 

Kielkoof &; Allard (1995^ 


1396.40 

2p<T„-3d(T g 

ves (Koester et al. 1985, e.g.) 

e.e. Allard et al. f 19991^1 

Lyfl 

1035.77 

lscr 5 -6h(7ii 

no 

Allard, Kielkoof & Feautrier ('1998'Pl 


1060.99 

lscr g -4f7T„ 

ves (Koester et al. 1996) 

Allard. Kielkoof & Feautrier 11998121 


1078.34 

2p<T u -5g7T g 

ves (Koester et al. 1996) 

Allard. Kielkoof &: Feautrier C19981S1 

Lyq 

978.24 

lso-„-8k(T„k 

no 

Allard, Noselidze &; Kruk (2009p^ 


992.05 

lscr 5 -8k(7 u 

ves (Wolff et al. 2001) 

Allard, Noselidze & Kruk (2009 


996.03 

2pcr u-7ln g 

ves (Wolff et al. 2001) 

Allard, Noselidze & Kruk (20091^1 

Ly<5 

953.03 

lscTg-lOmcTu 

2p<7 u -9l7T 5 

no 

Santos &; Kepler (2012) c 


961.61 

lscr 9 -8k7T u 

2p(7 u -91cTc, 

no 

Santos &; Kepler (20121^ 


a via unified theory 

^ The level 8k<7t4 is often miswritten in the literature as 8jcr u - The letter j, by convention, must not used 
for orbital representation. 
c via quasi-static approach 



Figure 2. Normalized profiles for the Lyman a line. In black 
for T = 20 000 K, dark grey for T = 15 000 K and light grey 
for T = 10 000 K. The dotted line shows calculations for per- 
turber density of 10 16 cm -3 , the continuous for 10 17 cm -3 and 
the dashed for 10 18 cm -3 . The same labels are adopted for the 
other Lyman lines. The smaller box shows in detail the behaviour 
of the satellite at ~ 1400 A with changes in temperature for den¬ 
sity of 10 17 cm -3 . 


data for the stars HS 0507+0434A and WD 0231-054 in the 
region of Lyman a and calculated models using the codes 
Tlus ty (jHnbenvl Il988h and Synspec llHubenv fc Land 
1201 j) . The dashed line shows the model obtained using the 
quasi-molecular contributions calculated in this work, and 


the dotted line uses tables provided by N. F. Allarcf). The 
models are not fitted to the data, but simply over-plotted 
considering normalization at an arbitrary wavelength and 
assuming effective temperature and surface gravity found in 
the literature. In both cases, there are two major differences 
between the models: the blue wing of the central line and the 
position of the satellite. Allard’s models seem to fit better 
the central line, but the central wavelength of the satellite 
is underestimated, being better determined by our model. A 
possible explanation for such difference is that Allard’s cal¬ 
culations use potential energy data only to small distances, 
so the real asymptotic value of each level is not reached. 
As a subtraction of the calculated asymptotic value [so that 
V(R max ) = 0] is necessary to the numerical evaluation of 
the line, a larger value than the real one is subtracted, over¬ 
estimating the shift in the difference potential and making 
the satellite appear closer to the line centre. Such problem 
does not happen in our calculations, since we evaluated the 
potential energy data to large distances. 


3.2 Lyman /3 

Ten transitions contribute to the Lyman /3 profile. Their 
combined effect gives the total profile shown in Fig. [5] Our 
calculations are shown for temperatures of 10 000, 15 000 
and 20 000 K and perturber densities of 10 16 , 10 17 and 
10 18 cm -3 . Apart from the strong satellites at 1060 and 


2 http://mygepi.obspm.fr/~allard/ 
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HS 0507+0434A 



1150 1200 1250 1300 1350 1400 

MA) 


Figure 3. HST spectrum of the white dwarf HS 0507+0434A 
(grey) and TluSTy/Synspec models using the quasi-molecular 
contribution tables evaluated by us (dashed line) and by N. 
F. Allard (dotted line), assuming temperature o f 20 000 K and 
log(g) = 7.75 [close to the values determined by [Andrews et ahl 
j2012T) l. In the small box, the spectrum was smoothed to enhance 
the position of the satellite. 



MA) 


Figure 4. HST spectrum of the white dwarf WD 0231-054 (grey) 
and Tlusty/Synspec models using the quasi-molecular contri¬ 
bution tables evaluated by us (dashed line) and by N. F. Allard 
(dotted line), assuming temperature of 15 000 K and log( g) = 8.80 
[simil ar to the values determined bv lGianninas, Bergeron &; Ruid 
d201lh |. Again we show the smoothed spectrum in the satellite 
region to enhance the satellite position, which is better fitted by 
our data. 


1078 A already studied in white dwarf spectra, there is an¬ 
other satellite around 1035 A, probably too close to the line 
centre to be identified in white dwarfs, since the broaden¬ 
ing of the line due to other mechanisms hides it. Temper¬ 
ature is again not an important factor, result that is cur¬ 
rent in the literature. However, again in a more detailed 
scale, the strong satellites do show some variation. They 
also show small oscillations. Possible reasons for that are 
the fact that the calculations presented here assume that 
the perturbers have a single constant velocity and no aver¬ 
aging was do ne, cau s ing a n effect similar to diffraction, as 
mentioned bv lAllar ara and, mainly, the contribution of 



MA) 


Figure 5. Normalized profiles for the Lyman (3 line. The smaller 
box shows in detail the behaviour of the two stronger satellites, 
at 1060 and 1078 A, with changes in temperature for density 
10 17 cm -3 . 



MA) 


Figure 6. Normalized profiles for the Lyman 7 line. The smaller 
box shows in detail the behaviour of the satellites at 990 and 
995 A with changes in temperature for density 10 17 cm -3 . 

the window fu nctio n in the Fou rier trans form [as reported 
bv 1 Allard. Sahal-Brechot fc Biraudl dl974h l. 

3.3 Lyman 7 

Fourteen transitions of the Hj quasi-molecule contribute to 
Lyman 7 . The calculated total profiles, for temperatures of 
10 000, 15 000 and 20 000 K and perturber densities of 10 16 , 
10 17 and 10 18 cnF 3 , are shown in Fig. [B] The noticeable 
satellite features are at 979, 990 and 995 A. The first one is 
again to close to the line centre for detection in stars, but 
the combined effect of the other two was already detected, as 
detailed in the introduction. Since this sate llit e has two com ¬ 
pon ents, as already reported by lAllard. Noselidze fe Krukl 
(1200911 . its behaviour with temperature shall not be very 
simple, because each component has its independent varia¬ 
tion. The rest of the line shows no significant change. 

In Figs. [7] and [5] we show FUSE observed spectra in 
the region of the lines Lyman /3 and Lyman 7 compared to 
calculated models using the codes Tlusty (lHubenvlll988li 
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MA) 


Figure 7. FUSE spectrum of the white dwarf Wolf 1346 (grey) 
and TluSTy/Synspec models using the quasi-molecular contri¬ 
bution tables evaluated by us (dashed line) and by N. F. Allard 
(dotted line), assumi ng temperature of 2 0 700 K and log(g) = 8.00 
[as determined bv lGiammichele, Bergeron, fc Dufourl i2012l ll. 



MA) 


Figure 8. FUSE spectrum of the white dwarf CD- 
38 0 10980 (grey) and TluSTy/Synspec models using the 
quasi-molecular contribution tables evaluated by us (dashed 
line) and by N. F. Allard (dotted line), assuming temper- 
ature of 26 000 K and l og(g ) = 8.00 [as determined by 
iGiammichele, Bergeron, fe Dnfoiirj li2012h [. 


and Synspec dHubenv fc LanJl201ll ). As before, the dashed 
line shows the model which uses the quasi-molecular con¬ 
tributions calculated in this work and the dotted line uses 
tables by N. F. Allard. In this case, there seems to be no 
difference between the central line of each model, but the 
position of the satellite is again underestimated by Allard’s 
models, probably due to the same reason as for the Lyman a 
satellite. The underestimation becomes more evident as one 
moves to high-order lines because higher levels converge at 
larger distances, so that the higher the line, the larger is 
the error in the calculated asymptotic value and, by con¬ 
sequence, the larger the shift of the satellite towards the 
central line. 



MA) 


Figure 9. Normalized profiles for the Lyman (5 line. The smaller 
box shows in detail the behaviour of the satellite at 962 A with 
changes in temperature for density 10 17 cm -3 . 


3.4 Lyman <5 

Lyman 8 has contributions of 18 transitions. It is evaluated 
with the unified theory for the first time in this paper. The 
normalized profiles are shown in Fig. [5] for temperatures of 
10 000, 15 000 and 20 000 K and perturber densities of 10 16 , 
10 17 and 10 18 cm -3 . We find stronger satellites around 953 
and 962 A. Neither has reported detection in the literature. 
The former probably is to close to the line centre for that, 
but the latter, although not mentioned in the literature, is 
present in the FUSE spectrum of CD -38° 10980 published 
in fig. 4 of IWolff et al.l (1200 il l . Similar to what was noted 
for the 1400 A satellite in Ly a, we see that its peak moves 
to smaller densities and that it broadens as temperature 
increases, so its behaviour is probably also a good diagnostic 
for temperature determinations. 


4 LINE PROFILES FOR THE BALMER 
SERIES 

The Balmer series is composed by transitions that have 
n = 2 as the final level, in case of absorption, or the ini¬ 
tial level, for emission. It is located in the visible portion of 
the spectrum, between 3646 A (H oo) and 6563 A (H a). We 
evaluate here every line from HatoH 10, as can be seen in 
the opacity profiles in Fig. 1101 for temperature of 10 000 K 
and perturber density of 10 17 cm -3 . The stronger satellites 
are listed in Table [2] Only one of them was observed so far, 
which is probably due to lack of laboratory data away from 
the line centre for the other lines. Such satellites would be 
stronger for temperatures around 10 000 K, but the ioniza¬ 
tion rate of hydrogen is low at such temperature, so they are 
effectively stronger around 15 000 K, when the H + popula¬ 
tion is more significant. In white dwarfs, lines above H 11 are 
hardly ever seen, since the high pressure causes then to over¬ 
lap and disappear. The profiles evaluated here are probably 
sufficient to study any spectra in the visible region. 
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Table 2. Wavelengths for the maximum of the stronger satellites (if existent) in each calculated line in the 
Balmer series (for T = 10 000 K and n = 10 17 cm -3 ), the transition(s) that originates them and notes on 
observation and previous theoretical models. 


Line 

Satellite (A) 

Transition (s) 

Observed? 

Previous model 

Ho 

8506.60 

3d7r 5 -4f7Tu 

4ia u -bgn g 

ves (Kielkopf, Allard &; Decrette 2002) 

Kielkoof. Allard & Decrette (2002) a 

H/3 

5413.33 

3d7r£-6h7ru 

no 

Santos &; Kepler (2012) ^ 

H 7 

4392.21 

3 d 7 r 9 - 8 k< 5 u 
3d?r 9 -9k a u 
3pcr u -7ga g 

no 

Santos & Kepler (2012^ 


4571.88 

2pa u -7\& g 

3 der g - 8 k 7 r u 

3d7r g -8k7r u 

no 

Santos & Kepler (2012'P^ 

H<5 

4203.87 

2 pcr u -101(Jc, 

3d7T 5 -9k7T u 

no 

Santos & Kepler (2012^ 

He 

4005.35 

2 s< 70 -lOk 7 T ti 
2 pcr u - 10 i 6 g 
2pcr u -lll Sg 
2pa u -12na g 
3pcr u -12n (T g 
3d7r ff -10k7r u 
4fcr u -12ncr g 

no 

Santos & Kepler ('2012'JS 


via unified theory 

via quasi-static approach 



Figure 10. Opacity profiles for every computed line in the 
Balmer series for a temperature of 10 000 K and perturber density 
of 10 17 cm -3 . The satellites are not very intense in this case, but 
it is worth noticing that all lines are very asymmetrical, especially 
the high order ones (in the detailed box), in which the red wing 
falls of very quickly, while the blue wing has a smother behaviour. 


4.1 Balmer a 

For the Balmer a line, there are contributions due to 32 tran¬ 
sitions. The total profile is shown in Fig. [IT] for a tempera¬ 
ture of 10 000 K and densities of 10 16 , 10 17 and 10 18 cm -3 . 


The strongest satellite is a round 8500 A and has b een ob¬ 
served in plasma spectra bv lKielkopf fe Allardl (Il995l ). How¬ 
ever, in white dwarfs, it has not been identified for two rea¬ 
sons: first, it is too close to the region were the Paschen 
series is and, second, because, for temperatures at which 
the population is relevant (10 000-25 000 K), the star 
emits predominantly in the blue, having little flux in the 
region of this satellite. The satellite also contributes to the 
line asymmetry. 

4.2 Balmer /3 

Forty-six transitions contribute to Balmer /3. Its total profile 
is shown in Fig. 1121 for a temperature of 10 000 K and den¬ 
sities of 10 16 , 10 17 and 10 18 cm -3 . We can see many faint 
satellites in the red wing, causing asymmetry. The strongest 
satellite is around 5400 A, but it is probably too faint com¬ 
pared to the line intensity to be observed in white dwarfs. 
It is important to notice, though, that it contributes to the 
asymmetry, which, is worth remembering, does not appear if 
the interaction between hydrogen and proton is treated sim¬ 
ply as Stark effect, as is done in current atmosphere models. 

4.3 Balmer 7 

To Balmer 7 there are contributions due to 60 transi¬ 
tions. The total profile is shown in Fig. [13] for a temper¬ 
ature of 10 000 K and perturber densities of 10 16 , 10 17 and 
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Figure 11. Normalized profiles for the Balmer a line. The tem¬ 
perature is 10 000 K and densities are 10 16 (light grey), 10 17 (dark 
grey) and 10 18 cm -3 (black). The same labels are adopted for 
the other Balmer lines. Various satellites appear in the red wing, 
making the line very asymmetric. 



Figure 12. Normalized profiles for the Balmer (3 line. The 
stronger satellite is at 5400 A, but it is still too faint for observa¬ 
tion, its main importance being thus the asymmetry it causes. 


10 18 cm 3 . There are two noticeable satellites, close to 4400 
and 4600 A. The former is too close to the line centre to 
be observed in white dwarfs, the latter, despite being away 
from the centre, has very small intensity over the contin¬ 
uum, what would difficult or even prevent its observation. 
However, both of them cause asymmetry. 



MA) 


Figure 13. Normalized profiles for the Balmer 7 line. There are 
two strong satellite features, at 4400 and 4600 A, none of them 
reported in observed spectra until now. 



Figure 14. Normalized profiles for the Balmer 8 line. The most 
noticeable satellite is at 4200 A, but again it seems too faint to 
be observed in white dwarf spectra. 


4.5 Balmer e 

There are 88 transitions that contribute to Balmer t. We 
show the total profile in Fig.[l5]for a temperature of 10 000 K 
and perturber densities of 10 16 , 10 17 and 10 18 cm -3 . We 
can see a very faint satellite around 4005 A, with intensity 
of the same order of the continuum in the blue wing of the 
line, meaning it is probably not observable. However, as the 
intensity of the red wing rapidly decreases after the satellite, 
while the intensity in the blue line changes slowly, knowing 
the satellite is important to model the line asymmetry. 


4.4 Balmer 8 

We have to account for contributions of 74 transitions in 
Balmer 5. Its profile for a temperature of 10 000 K and per¬ 
turber densities of 10 16 , 10 17 and 10 18 cm -3 can be seen in 
Fig- El The strongest satellite is close to 4200 A, but, as for 
the 4600 A satellite in Balmer 7 , its intensity does not seem 
to be strong enough, when compared to the continuum, to 
allow its detection in white dwarf spectra. 


4.6 Balmer 8 , Balmer 9 and Balmer 10 

The lines of Balmer 8 , Balmer 9 and Balmer 10 have con¬ 
tributions due to 102, 116 and 130 transitions, respectively. 
Their profiles are shown in Fig. 1161 Fig. [T7] and Fig. 1181 
again for a temperature of 10 000 K and perturber densities 
of 10 16 , 10 17 and 10 18 cm -3 . It is worth noticing that, for 
n = 10 18 cm -3 , these high order lines start to fade due to the 
high pressure. There are no satellites visible in those lines. 
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Figure 15. Normalized profiles for the Balmer e line. A satellite 
feature appears at 4005 A, and the intensity of the profile drops 
after it. 



Figure 17. Normalized profiles for the Balmer 9 line. In spite of 
the absence of satellites, we can note a line asymmetry that does 
not exist in pure-Stark profiles. 



Figure 16. Normalized profiles for the Balmer 8 line. No satellites 
are seem, but the line is largely asymmetric. 

One might wonder in that case if the quasi-molecular con¬ 
tributions really matter. However, it is very clear that they 
do affect the line behaviour, again causing a strong asym¬ 
metry between the red and blue wings, which is until now 
completely ignored in white dwarf stars atmosphere models. 


5 COMPARISON TO ELECTRONIC STARK 

In current synthetic spectra for DA white dwarf stars, the 
effect of the interaction between free protons and radiating 
hydrogen is assumed to be exactly the same of the inter¬ 
action of free electrons with hydrogen. So the profiles are 
calculated with density n e of electrons, and such density 
is simply doubled to account for protons as well, as the 
number of free protons must be the same as the number 
of free electrons, because both originate, in this case, when 
hydrogen is ionized. As already said in the introduction, such 
approach disregards the quasi-molecular contributions that 
appear when temperature and pressure are adequate to the 
formation of a transient bond state between hydrogen and 
proton. 

This is not a good practice and must be affecting the 



Figure 18. Normalized profiles for the Balmer 10 line. Again no 
satellites, but a large asymmetry between red and blue wings can 
be seen. 

determination of temperature and surface gravity. We plot 
here the calculated lines for T = 10 000 K and log (n) = 17 
together with the correspondent electronic Stark effect to 
show that. Profiles for the Lyman series are shown in Fig. 
m and for the Balmer series in Fig. [20] It is easy to no¬ 
tice two large differences: the existence of satellites and the 
large asymmetries between red and blue wings. The former 
is more important for low order lines (n/ ^ 5), when the 
red wing usually shows at least one strong satellite. The lat¬ 
ter is evident for high order lines (ny > 5), in which the 
red wing falls quickly for the quasi-molecular profile, while 
the red wing in the electronic Stark behaves similarly to the 
blue wing and decreases slowly. These differenc es may b£_the 
cause of th e inconsistent densities report e d by Falcon et al.l 
(|2014h and iGenest-Beaulieu fe Bergeronl (12014 1. which are 
obtained when different lines are fitted. 

It is worth mentioning that the validity domain of the 
VCS theory, in which the Stark model s sho wn here are 
based, is Auj/uj p ^ 1 (IVidal, Cooper fe Smithlll97llf . where 
Alo is the angular frequency measured from the line cen¬ 
tre and uj p = 47rn e e 2 / m e is the plasma frequency (cgs 
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Figure 19. Comparison between each calculated line in the Ly¬ 
man series ( black continuous line) and the correspondent pure 
electronic Stark profile ( grey dashed line) for T = 10 000 K and 
log(n e ) = log(n p ) = 17. There is no great difference if one consid¬ 
ers only the blue wing, but the red wing is very different, show¬ 
ing satellite features when the quasi-molecular contributions are 
taken into account. 




3800 3900 4000 4100 



Figure 20. Comparison between calculated lines in the Balmer 
series ( black continuous line) and the correspondent pure elec¬ 
tronic Stark profile ( grey dashed line) for T = 10000 K and 
log(n e ) = log(rip) = 17. For the lower order lines, we see that the 
satellite in the blue wing must have a considerable contribution. 
On the other hand, in the high order lines, the major difference is 
the rapid decreasing of the red wing compared to the blue wing, 
effect that does not appear in the electronic Stark profiles. 


units). According to that, the profiles would be valid only in 
the central 1-10 A of the line. However, the validity of the 
theory has been extrapolated in the current models, since 
observed lines of white dwarf stars are usually more than 
100 A wide. 


is only appropriate for the frequencies in the line wings. 
Therefore, the central b ehaviour of the profiles evaluated by 
ISantos fe Keplei] (j2012l ) should not be relied upon. Besides 
that, such profiles were calculated summing up contributions 
due to different transitions, when, actually, they should be 
convoluted, since lifetimes for hydrogen excited states are of 
the order of 10 -8 s, while the typical duration of collisions 
in the physical conditions of white dwarfs is around 10 -13 - 
1 CT 12 s, implying that different states may form throughout 
the electronic transition and contribute to the line. 

The profiles calculated here are based on the unified 
theory, in which no assumption is made concerning the du¬ 
ration of collisions, and the partial profiles are convoluted, 
so the form of the line can be relied upon for the whole wave¬ 
length range. Nevertheless, we recall that both the unified 
theory and the quasi-static approach assume adiabatic per¬ 
turbations, disregarding interference between different lev¬ 
els ,_which_can_result in underestimation of the line shift 
dAllard et al.lll999h . However, as the behaviour of the line 
centre in white dwarf stars is dominated by electronic Stark 
and Doppler effects, being the quasi-molecular contributions 
more important in the wings, the profiles shown in the pre¬ 
vious sessions should be satisfactory for model atmosphere 
applications. 

Ho wever, wit h the adve nt of experiments such as the one 
done bv lFalcon et al.1 (l2012l l at the Z Pulsed Power Facility 
at Sandia National Laboratories, more accurate profiles in 
the central region of the line are necessary. When trying to 
evaluate such de t ailed profiles using the published data of 
ISantos fc Keplerl (l 2012 l h we noticed that the data was not 
accurate enough for that. There were discontinuities in some 
of the potential data, due to convergence issues in the im¬ 
plemented computational methods, as can be seen on Fig. 
EU That lead to abnormal behaviour in the line centre, as 
shown for Lyman a in Fig. [22] To avoid that behaviour, it 
was necessary to truncate the potentials at ~180 ao, where 
ao « 0.5292 A is the Bohr radius, but that way the asymp¬ 
totic behaviour is not reached and the line centre is still not 
accurately determined (Fig. 1231) . since it is formed mainly 
when the particles are far apart. Besides, as the mean dis¬ 
tance between the particles for density of 10 17 cm -3 , for ex¬ 
ample, is around 250 ao, going up to 180 ao is not sufficient 
if one is concerned with the line centre. As at that distance 
the asymptotic behaviour of the potential energy, which is 
oc r -2 , is still not reached, it is also not possible to extrapo¬ 
late the data with the necessary precision. For larger values 
of n, the asymptotic behaviour is reached even farther away, 
so the truncation may largely affect the determination of the 
line centre profile. 

Thus, we recalculated the potential data for the H^ 
quasi-molecule. In order to check the convergence of the po¬ 
tentials to their asymptotic values, we had to make use of 
a larger number of terms in the recurrence relation^], what 
leads to a larger set of homogeneous linear equations and 
demands, in consequence, larger matrices to solve. We cal- 


6 THE LINE CENTRE 

Quasi-static profiles are inherently a poor approximation to 
the line centre, since its main assumption of zero velocity 
implies that only collisions whose duration is larger than 
the interval between them are taken into account, and that 


3 In ISantos fc Kepler] (1 20 1 2j ). the recurrence relations were trun¬ 
cated at indexes equal to 25. Now, for n=l and n=2 for example, 
we had to truncate the indexes at 300 to get the correct asymp¬ 
totic value to a precision at the 21 st decimal. 
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Figure 21. Difference potentials (minus the asymptotic value) for 
the six transitions contributing to the Lyman-a line. We can note 
that in the old data there is a discontinuity, which, transformed 
to wavelength, corresponds to >1.2 A. This issue was solved in 
the new data. 



Figure 22. Profile for the centre of the Lyman a line for a tem¬ 
perature of 10 000 K and pertur ber density of 10 1 ' cm~ 3 evalu¬ 
ated using the data published in lSantos fe Kepler (2012) (dashed 
line) and the profile calculated with the improved data ( continu¬ 
ous line). We can note an abnormal behaviour in the old profiles, 
which is a consequence of discontinuities in the potential data. 



MA) 


Figure 23. Same as Fig. 1221 but now truncating the old po¬ 
tentials before the discontinuity. The abnormal behaviour disap¬ 
pears, but as the asymptotic behaviour of the potential is not 
reached, the contributions to the centre are not so well deter¬ 
mined. A better determination is obtained with the improved 
data. 



Figure 24. Profile for the centre of the Lyman /3 line for a tem¬ 
perature of 10 000 K and pertur ber density of 1 0 17 cm" 3 evalu¬ 
ated using the data published in lSantos fe Kepler) i2012l ) (dashed 
line) and the profile calculated with the improved data ( continu¬ 
ous line). The broadening of the line is underestimated when one 
uses the old data, affecting also the position of the satellite. 


culated the profiles for Lyman a, Lyman /3 and Balmer a 
with the new data. 

The obtained profile for Lyman a is shown in both Fig. 
[22]and Fig.[23]to comparison with the profiles with previous 
data. It is extremely different than the one calculated with 
the old data without truncation, since it has no discontinu¬ 
ities (Fig. 1221) . Compared to the profile calculated with the 
old data truncated in 180 ao (Fig. 12311 . the difference is little, 
but it still important, especially in the central wavelengths, 
were we can see that the broadening is underestimated if 
one does not consider the behaviour at large distances. The 
same result is obtained when one analyses the profiles for 
Lyman /? (Fig. 1241) and Balmer a (Fig. 1251) . with the ad¬ 
diction that it becomes apparent that such underestimation 
may also affect the position of the satellite, which will ap¬ 
pear closer to the line centre. 


7 CONCLUSIONS 


The quasi-molecular contributions of Hj are an effect very 
important to the shape of the lines in a hydrogen gas. The 
most remarkable contributions are the satellite lines, which 
are a consequence of the formation of transient bond states 
between hydrogen atoms and free protons due to an ex¬ 
tremum in the potential difference. However, that is not 
the only effect of the quasi-molecular absorptions: they are 
also responsible for large asymmetries in the line wings. 
The satellit es are an effec t often discussed in the litera¬ 
ture (e^^ja^^^^oestei] J 


19921 : lAllard. Noselidze fe Krukl 


120091 : ISantos fe Kepler 2012), since they are quite easily ob¬ 


served in the ultraviolet. Meanwhile, the asymmetry seems 
to be ignored, as it may not be directly detected. However, 
the fact that the quasi-molecular profile red wing falls down 
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Figure 25. Profile for the centre of the Balmer a line for a tem¬ 
perature of 10 000 K and pertur ber density of 1 0 17 cm^ 3 evalu¬ 
ated using the data published in lSantos fc Kepler! ||2012| ) ( dashed 
line) and the profile calculated with the improved data ( continu¬ 
ous line). Again it is clear that the broadening is underestimated 
when the potential is truncated. 

very quickly, especially for transitions involving states with 
high n, is probably very important in the opacity determi¬ 
nations. Not taking it into account, and continue to use the 
symmetric Stark profiles for protons, is almost certainly af¬ 
fecting not only the determinations of temperature, but also 
of gravity, to which such high order lines are more sensible. 

The only line considered at present in this sense is Ly¬ 
man a, whose quasi-molecular contributions turned out to 
be very important to the opaci ty of cool white dwarfs (T e g < 
6 000 K), as shown firstly by (Kowalski fe Saumonl (1200611 . 
It is worth noticing that, at such temperatures, satellite 
lines are not even observed, but the quasi-molecular effect 
is important anyhow. The fir st step towards g e nerali zation 
of those results was done by ISantos fe Kepler! (l2012h . who 
calculated all potentials energies for n up to 10 and dipole 
moments for the Lyman series up to Ly 8 and the Balmer 
series up to H 10. They also evaluated the quasi-static pro¬ 
files for those lines, and already pointed out there that the 
quasi-molecular effects in the line wings shall probably be 
important for opacity, but remain unaccounted for in the 
current models, which assume the interaction between hy¬ 
drogen and protons is completely analogous to that of hydro¬ 
gen with electrons. Considering that the quasi-static profiles 
are not correct at the line centres, we re-evaluated here all 
those line profiles, now using the unified theory, that can give 
the whole line shape. The data files are available to the com¬ 
munity as online material accompanying this paper and also 
at our web page (http://astro.if.ufrgs.br/ingrid/index.htm). 
We believe that such profiles are enough to study both Ly¬ 
man and Balmer series in white dwarfs. 

We also noted that the precision of the current poten¬ 
tial energy data was not enough to compare the calculated 
profiles to experimental results. The main reason for that 
were discontinuities due to numerical imprecision. We solve 
that in the code by increasing the truncation in our recur¬ 
rence relations from index 25 to 300. Results were shown 
in section [S] for Lyman a, Lyman /3 and Balmer a and a 
large improvement could be verified. We hope to be able to 
compare our profiles to experimental data and actually con¬ 


firm the importance of quasi-molecular contributions due to 
H^" even when no satellites are observed, in order to make 
their inclusion in atmosphere models more consistent and 
consolidated. 
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